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ABSTRACT 

We present the results of high-resolution AP'^M+SPH simulations of merging clus- 
ters of galaxies. We find that the compression and shocking of the core gas during a 
merger can lead to large increases in bolometric X-ray luminosities and emission- 
weighted temperatures of clusters. Cooling flows are completely disrupted during 
equal-mass mergers, with the mass deposition rate dropping to zero as the cores of the 
clusters collide. The large increase in the cooling time of the core gas strongly suggests 
that cooling flows will not recover from such a merger within a Hubble time. Mergers 
with subclumps having 1/8'^ of the mass of the main cluster are also found to disrupt 
a cooling flow if the merger is head-on. However, in this case the entropy injected into 
the core gas is rapidly radiated away and the cooling flow restarts within a few Gyr 
of the merger. Mergers in which the subcluster has an impact parameter of 500kpc do 
not disrupt the cooling flow, although the mass deposition rate is reduced by ~ 30%. 
Finally, we find that equal mass, off-centre mergers can effectively mix gas in the cores 
of clusters, while head on mergers lead to very little mixing. Gas stripped from the 
outer layers of subclumps results in parts of the outer layers of the main cluster being 
well mixed, although they have little effect on the gas in the core of the cluster. None 
of the mergers examined here resulted in the ICM being well mixed globally. 
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1 INTRODUCTION 

In the Cold Dark Matter scenario, clusters of galaxies form 
through a succession of mergers of smaller subclusters and 
groups. Major mergers can be the most energetic events 
in the Universe since the Big Bang, with hydrodynamic 
shocks dissipating much of the kinetic energy of the col- 
lision (10''^"^"' ergs) into the Intracluster Medium (ICM). 
The theoretical picture is supported by X-ray observations, 
which suggest that many clusters of galaxies are either un- 
dergoing mergers or have experienced merger events in the 
recent past. ROSAT observations generally focussed on dis- 
tortion of the X-ray isophotes (e.g. White, Briel & Henry 
1993) and the presence of X-ray substructure (Forman & 
Jones 1994; Buote & Tsai 1996), while evidence for shocks 
in the ICM was seen in temperature maps created from 
ASCA observations (e.g. Arnaud et al. 1994). The launch of 
Chandra and XMM-Newton has provided a wealth of obser- 
vational data, including high-resolution temperature maps 
of shocks in merging clusters (Markevich & Vikhlinin 2001) 
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and the surprising discovery of sharp surface brightness dis- 
continuities between the hot cluster gas and moving "cold 
fronts" (Vikhlinin, Markevitch & Murray 2001), which are 
thought to be the cores of subclusters which have survived 
the merger process (Markevitch et al. 2000). 

Modern cosmological simulations (Bertschinger 1998 
and references therein; see also Frenk et al. 1999) are ca- 
pable of following the evolution of large volumes of the uni- 
verse from early times until the present day, but are not 
particularly suited to modeling mergers as the resolution of 
individual clusters is limited unless techniques such as res- 
imulation (e.g. Eke, Navarro & Frenk 1998) are used. Cos- 
mological simulations are also hard to analyse as the clus- 
ter does not necessarily have time to return to hydrostatic 
equilibrium between mergers, making it hard to separate 
one merger from the next. An alternative is to study the 
merger of idealized systems, which, while less realistic than 
the full cosmological approach, have a number of advan- 
tages. Firstly, the objects participating in the merger are 
fully formed, and so there is no need to simulate the large 
volumes of space necessary to capture the evolution of a 
cluster of galaxies from high redshift. We can therefore use 
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a smaller volume in which the merging systems are well re- 
solved, and this also has implications for our hydrodynamic 
scheme; this is discussed in Section H. As we are studying 
a single merger event between two objects which are ini- 
tially in hydrostatic equilibrium changes in the final object 
are easily quantified, and the controlled initial conditions al- 
low parameter space to be explored. While full cosmological 
simulations are ultimately more realistic, the simplified pic- 
ture represents the best chance of understanding the physics 
underlying the merger process. 

Several authors have presented results from such sim- 
ulations. Head-on mergers between King-model clusters of 
unequal mass containing gas and dark matter were studied 
by Roettiger, Burns & Loken (1993) and Roettiger, Loken & 
Burns (1997, hereafter RLB97) using a finite-difference code 
to model the hydrodynamics and an A''-body code to model 
the evolution of the gravitational potential, while Schindler 
& Miiller (1993) used a similar approach, although with 
less controlled initial conditions. However, these simulations 
treated the baryonic component as being massless, limit- 
ing the results to very gas-poor clusters. Pearce, Thomas 
& Couchman (1994, hereafter PTC94) examined head-on 
mergers of systems containing gas and dark matter using an 
adaptive particle-particle, particle-mesh (AP^M; Couchman 
1991) code with the gas being modeled using Smoothed- 
Particle Hydrodynamics (SPH; Monaghan 1992), with their 
results focusing mainly on the evolution of the core of the 
cluster rather than the observable properties of the merger. 

More recently Roettiger, Stone & Mushotzky (1998) 
have used a hydrodynamic code based on the Piecewise- 
Parabolic Method (PPM; CoUela & Woodward 1984) cou- 
pled to a Particle-Mesh A'^-body code to examine the evolu- 
tion of the cluster A754, suggesting that its unusual X-ray 
morphology (e.g. Henry & Briel 1995) is the result of a recent 
(< 0.5Gyr), slightly off-centre merger. Ricker (1998) also 
used PPM to simulate off-centre cluster mergers, although 
these simulations do not include the dark matter which typ- 
ically dominates the dynamics of clusters and can make a 
significant difference to the end state of the merger pro- 
cess (PTC94; see also Section |!|). Finally, Takizawa (1999; 
2000) carried out simulations of merging clusters using SPH, 
with gravitational forces calculated using a tree algorithm 
(Barnes & Hut 1986). The simulations of PTC94 and Tak- 
izawa (1999) are the most similar to our own and provide a 
useful comparison, although both were carried out at lower 
resolution than the results described here. 

In this paper we present results from high-resolution 
(131072 particle) simulations of the merger of idealized clus- 
ters of galaxies containing both dark matter and gas, using 
an AP'^M+SPH code. Off-centre mergers and mergers be- 
tween unequal-mass systems are also examined, with ref- 
erence to the observable properties of clusters of galaxies. 
The layout of the paper is as follows; in Section we de- 
scribe our hydrodynamic scheme and the generation_of ini- 
tial conditions, while results are described in Section 
significance of the results are discussed in Section 
conclusions are drawn in Section H. 




2 METHOD 

The simulations described here have been carried out using 
HYDRAQ (Couchman, Thomas & Pearce 1995), an A-body, 
hydrodynamics code which combines an AP^M A-body al- 
gorithm (Couchman 1991) with SPH. Tests of our code can 
be found in Couchman et al. (1995), Thacker et al. (2000) 
and Ritchie & Thomas (2001). 

In SPH, the equations of motion for a compressible fluid 
(e.g. Landau & Lifshitz 1959) are solved using a Lagrangian 
formulation in which the fluid is partitioned into elements, a 
subset of which are represented by particles of known mass 
m and specific energy en. Continuous fields are represented 
by interpolating between particles using a smoothing kernel, 
which is normally defined in terms of a sphere containing a 
fixed number of neighbours, centred on the particle in ques- 
tion. The radius of the smoothing sphere is adjusted so as to 
keep the neighbour count approximately constant, making 
SPH adaptive in both space and time. The particle nature of 
SPH means that there is no grid to constrain the geometry or 
dynamic range of the system being studied, and allows SPH 
to be easily integrated with many A-body solvers. However, 
unlike PPM, SPH requires an artificial viscosity to convert 
relative motion to heat; we use a pairwise artificial viscos- 
ity (Monaghan & Gingold 1983) as described by Thacker et 
al.(2000). 

SPH is often criticised for capturing shocks poorly in 
comparison with modern high-order Godunov-type schemes, 
and for having poor resolution in low-density regions, which 
are represented by relatively few particles. While both of 
these criticisms can be true, they are of limited importance 
for the simulations considered here. Steinmetz & Miiller 
(1993) find that the shock capturing ability of SPH is closely 
linked to the number of particles in the system being mod- 
eled, with SPH giving accurate results if in excess of 10* par- 
ticles are used to model three-dimensional problems. This 
figure is out of reach of most cosmological simulations, but 
is easily achievable with our simplified approach, and the 
initial conditions are generated with the constraints of Stein- 
metz & Miiller (1993) in mind. We will not be able to follow 
the propagation of shocks into the lower-density outer re- 
gions of the cluster as well as finite-difference methods, but 
the X-ray luminosity of such regions will be low and so the 
limitations of our hydrodynamical method will not impact 
any potentially observable properties of the system. 

Radiative cooling is implemented in our code by adding 
a sink term ^ to the Energy equation, where ^ is the emis- 
sivity (the emission rate per unit volume), interpolated from 
the cooling function of Sutherland & Dopita (1993). Cool- 
ing during shock-heating (see Hutchings & Thomas 2000) 
is minimized by applying the artificial viscosity prior to the 
gas being allowed to cool. The cooling is assumed to occur at 
constant density, with the time-step ensuring that this con- 
dition is approximately satisfied, as described in Thomas 
& Couchman (1992). Particles which have cooled to lO'^K, 



^ This code is in the public d omain and can be downloade d from 
http://hydra.susscx.ac.uk/ or http://hydra.mcmaster.ca/ 



^ The specific energy t is related to the gas temperature T by 
e = 3kBT/2fimii, where fce is Boltzmann's constant, mn is the 
mass of the hydrogen atom and )i = 0.6 is the relative molecular 
mass. 
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Run no. 


Np 


Mdm(M0) 


Mgas{MxQ) 


po{cm ^) 


M1/M2 b/rcoTc 


ALx/Lx 


^Jew/Jew 


AE/E 


El 


131072 


1.25 X IQio 


1.13 X 


109 


10-3 


1 


0. 


7.99 


3.31 


5 X 10-* 


E2 


131072 


1.25 X IQif 


1.13 X 


10« 


10-3 


1 


2.5 


7.01 


3.17 


6 X 10-* 


E3 


131072 


1.25 X IQi" 


1.13 X 


lO'J 


10-3 


1 


5. 


5.01 


2.80 


1 X 10-* 


E4 


131072 


1.25 X IQi" 


1.13 X 


lO'J 


10-3 


1 


10. 


2.59 


1.86 


3 X 10-* 


Ul 


73728 


1.25 X 10"' 


1.13 X 


109 


10-3 


8 


0. 


2.26 


1.71 


8 X 10-* 


U2 


73728 


1.25 X IQi" 


1.13 X 


10« 


10-3 


8 


2.5 


1.95 


1.55 


9 X 10-* 


U3 


73728 


1.25 X IQio 


1.13 X 


109 


10-3 


8 


5. 


1.52 


1.32 


1.1 X 10-3 


CI 


131072 


1.25 X IQi" 


1.13 X 


109 


10-2 


1 


0. 


4.94 


4.34 


8 X 10-* 


C2 


73728 


1.25 X IQio 


1.13 X 


109 


10-2 


8 


0. 


2.18 


1.83 


4.4 X 10-3 


C3 


73728 


1.25 X IQio 


1.13 X 


109 


10-2 


8 


2.5 


1.09 


1.08 


3.9 X 10-3 


Rl 


32768 


5 X lO^o 


4.5 X 


109 


10-3 


1 


0. 


5.01 


3.34 


2.0 X 10-3 


R2 


8192 


2 X 10" 


1.8 X 


LO"' 


10-3 


1 


0. 


3.60 


2.94 


2.6 X 10-3 



Table 1. Details of the simulations examined in Section pi Listed are the run number, the total number of particles in the simulation 
(Dark matter + Gas), the mass of each dark matter particle, the mass of each gas particle, the core density, the ratio of masses of the 
two clusters M1/M2, the impact parameter b, the ratio of the peak bolometric X-ray luminosity to the bolometric X-ray luminosity 
6Gyr before the merger ALx/^x = Lt—o/Lt—_QQyr, the ratio of the peak emission-weighted temperature to the emission-weighted 
temperature 6Gyr times before the merger ATew/Tew = Tt=o/^t=— 6Gyr ^^nd the energy conservation during the simulation AE/E. 



the temperature at which the cooling function drops to zero, 
are converted to collisionless 'star' particles to avoid a build- 
up of cold gas that can cause the SPH algorithm problems 
(Pearce et al. 2001; Ritchie & Thomas 2001). We neglect the 
effect of thermal conductivity, which is known to be strongly 
suppressed in clusters of galaxies (e.g. Fabian, Nulsen & 
Canizares 1991; Vikhlinin et al. 2001). 



2.1 Initial Conditions 

Our idealized clusters initially have gas and dark matter 
density profiles given by a Hubble profile 



p{r) 



Po 



[1 + (r/rc)2]='''' ' 



(1) 



where po is the central mass density and Vc is the cluster 
core radius. This is similar to the gas density profile deduced 
from the X-ray surface brightness of clusters of galaxies, and 
has been widely used by other authors simulating cluster 
mergers (e.g. PTC94; RLB97; Kicker 1998; Takizawa 1999). 
A'^-body simulations (e.g. Moore et al. 1999) suggest that the 
dark matter is more centrally concentrated than that given 
by Equation hi, and Ricker & Sarazin (2001) have recently 
carried out simulations in which the dark matter density 
follows the Navarro, Frenk & White (1997) profile 



p{r) = 



Ps 



(r/r,){l+r/rs)^' 



(2) 



where ps and r^ are a scaling density and radius that are 
dependent on the halo mass. They find that in their simu- 
lations the morphological changes, temperature jumps and 
gas velocities that are similar to previous work in which the 
dark matter density is given by Equation |l|, although the 
peak X-ray luminosity is higher. 

In our simulations, we truncate the density profile at 
R = 16rc, a numerical compromise which allows the core to 
be well resolved yet be located well within the outer bound- 
ary of the system. To set up the density profile, we first place 
gas particles randomly within a cubical simulation volume. 



which is then evolved at constant temperature until spuri- 
ous fiuctuations arising from the initial particle distribution 
have died away. Particles are then ordered in terms of their 
distance from the centre of the box and are translated ra- 
dially to match the desired mass profile, ensuring that the 
gas starts off close to a relaxed state. Finally, collisionless 
dark matter particles with 9 times the gas mass are placed 
on top of each gas particle. 

Once particles have been placed, particle velocities must 
be set. Dark matter velocities are drawn from a 3-D Gaus- 
sian distribution of width a in each direction, where a is 
determined by solving the Jeans equation 



1 d(pa^) 
p dr 



GM 



(3) 



Gas particles have their initial velocities set to zero and are 
given a temperature equivalent to a^, so that 



/? 



1, 



(4) 



kBT/iirriH 

with particle densities set from Equation hi The sound cross- 
ing time for the cluster is 



Ca 



'yksTc 



= 1.75Gyr. 



(5) 



Figure hj shows the evolution of the density profile over 
a period of 20 sound crossing times for an isolated system 
set up with the density profile specified by Equation hi. Some 
expansion in the outer region due to the initially-truncated 
density profile is visible, but these fiuctuations are small and 
settle down within a couple of sound crossings. Our initial 
conditions for simulations of equal mass mergers simply take 
two such spheres, displaced by a distance equal to the impact 
parameter b in the y-direction and a distance 2R in the x- 
direction. The initial relative velocity in the a::-direction is 
equal to the circular speed 



^coll 



GM 
R ' 



(6) 



B. W. Ritchie and P. A. Thomas 



at the edge of the cluster, where M is the mass of each cluster 
and R = 16rc as above. This corresponds to a parabolic orbit 
but, once energy is dissipated in the encounter, leads to a 
bound system. The effect of varying the encounter velocity 
is discussed in PTC94. 

Mergers between systems of different mass are set up by 
using fewer particles in the lower-mass cluster (i.e. the mass 
per particle is kept constant) and scaling the second cluster 
to match the desired radius and velocity dispersion profile. 
We use the same impact parameter and relative velocity as 
in the equal mass encounters. 



2.2 Simulations 

Details of the simulations are listed in Table hi. With the ex- 



ception of the simulations carried out in Section 3.6, which 
examine the effect immerical resolution has on our results, 
all simulations of equal-mass mergers presented here use 
32768 particles each of gas and dark matter in each clus- 
ter (i.e. an equal-mass merger simulation contains a total of 
131072 particles). In simulations of unequal-mass mergers 
the particle mass is kept constant, so that the subcluster 
contains 1/A^ fewer particles, where A'^ = M1/M2 is the ra- 
tio of cluster masses. We are not attempting to simulate the 
evolution of clusters in a cosmological context, and therefore 
our simulations use a non-expanding volume with vacuum 
boundary conditions. Particles leaving the simulation vol- 
ume are removed, although the size of the box is such that 
very few particles escape by the end of the simulation. We 
use a gravitational softening r^ = 0.2rc, which is above the 
particle separation in the core. This choice of softening is 
discussed in detail in Pearce, Thomas & Couchman (1993). 
We use a gas particle mass of 1.25 x lO^M© and a dark 
matter mass 1.13 x lO^'^Mo, giving a total cluster mass of 
4 X IO'^'^Mq and virial temperature ~ 3keV. The subclus- 
ters used in simulating unequal mass mergers have a mass 
of 5 X IO^^Mq and virial temperature ~ 0.9keV. We take 
the core radius r^ to be lOOkpc, giving a central gas den- 
sity po ~ lO'^^cm""^, similar to that used by RLB97. The 
bremsstrahlung cooling time 



f'cool 



3 fceT 

2 MX) 



»^'>'"'°(TlK)'"(iF^)"-m 



(Sarazin 1986), where A(r) = AoiT/KY^'^ and Ao = 5.2 x 
10~^* erg cm"^ s~^, is therefore in excess of a Hubble time in 
the core of the cluster. While this might suggest that radia- 
tive cooling will have little impact on the results of our sim- 
ulations, cooling times can drop below a Hubble time when 
the density of the core gas increases as the cores collide and 
radiative cooling makes a significant (20 — 40%) dif fere nce 
to our results; this is discussed further in Section 3.3. In 



addition, we perform some simulati ons with a higher core 
density po ~ 10~'^cm~'^ (see Section 3.4) in which the cool- 
ing time is much less than a Hubble time, although icooi 
remains greater than the sound crossing time tac in all of 
our simulations. We therefore apply radiative cooling to all 
the simulations performed here. 

All of the simulations described here were carried out on 
a 700MHz Intel Pentium-3 workstation. The simulations of 
mergers between equal mass, low core density clusters take 
around 1800 timesteps and 15 hours of CPU time, while 
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Figure 1. Radial density profiles in a single-cluster test run. The 
initial density profile is plotted using a solid line, while the density 
profile after 20 sound crossing times is plotted using a dotted line. 
Both profiles have been normalised to po, the initial core density. 



simulations of mergers between high core density clusters 
take around 2500 timesteps and 21 hours of CPU time. 



3 RESULTS 

3.1 Morphology 

Figure ti shows a sequence of snapshots of the emission- 
weighted gas temperature (colours). X-ray surface bright- 
ness (contours) and velocity fields (arrows) during the head 
on merger between two systems with a masses of 4 x 10^*Mq 
(run El). The merging sequence is similar to that seen in 
other simulations (e.g. Ricker 1998; Takizawa 2000), and is 
therefore only described briefiy here. 

As the outer layers of the two clusters start to interact a 
weak shock forms (panels 1 & 2) , and gas is driven outwards 
in the plane perpendicular to the collision axis with the gas 
velocity increasing as it accelerates down the pressure gra- 
dient. As the cluster cores approach this shock strengthens 
(panel 3), as does the outfiow of gas in the plane of the 
collision. The increase in gas temperature between the two 
cores has become clearly visible in temperature maps of the 
system, but as the bulk of the X-ray emission is coming 
from the cores of the two clusters the shock still makes only 
a minor contribution to the integrated emi ssion- weighted 
temperature of the system (see Section 3.3). As the clus- 



ter cores collide (just before panel 4), a strong arc-shaped 
shock is driven into the outer layers of the cluster (visible 
in panel 4, and more clearly in panel 5). Meanwhile, gas 
in the core of the merged cluster goes through a period of 
expansion driven by the dark matter, cooling adiabatically 
to slightly less than the pre-coUision temperature with the 
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Figure 4. The final state of run E3, the 6 = 5 ofi^-centre merger 
shown in Figure pi Greyscales. contours and arrows represent the 
same quantities as in Figures ti and H, but have been rescaled to 
better show the rotation of the merger remnant. 



X-ray isophotes forming a bar shape along the collision axis, 
as noted by RLB97. Finally, the dark matter turns around 
and recoUapses to form a spherically symmetric final object 
(panel 6). Material ejected in the earlier stages of the merger 
can be seen falling back towards the core of the cluster. 

An off-centre merger with impact parameter 6 = 5 (run 
E3) is shown in Figure H. As before, the clusters have equal 
masses, and the behaviour is similar to that seen in Figure eI 
The shock that forms as the outer layers of the clusters in- 
teract is now oblique, but the shock and outflow is otherwise 
similar to the head-on merger. The picture changes some- 
what as the cluster cores interact, with the merger shock 
generated as the cores collide now propagating with a spi- 
ral pattern, and the cores partially surviving their first ap- 
proach (clearly visible in panel 5), completing most of an 
orbit before recoUapsing and merging completely. The final 
state of the cluster looks similar to that in the head-on case, 
although there is now significant rotation of the core, shown 
in more detail in Figure M. 



3.2 The structure of the merger remnant 

Figure H shows the radial gas density profiles of the merger 
remnant for three mergers with impact parameters 6 = 0, 2.5 
and 5 (runs El — E3) ten sound crossing times after the 
merger, when the gas has returned to hydrostatic equilib- 
rium and the profile has stopped evolving. The core den- 
sity is normalized to po, and the initial density profile is 
also plotted for reference. Best-fit values of the core radius 
Tc and core density pc are listed in Table y. The merging 
process has the effect of increasing the core radius and de- 
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Figure 5. Radial gas density profiles of the merger remnant 10 
sound crossing times after the merger, for three equal mass merg- 
ers with impact parameter b =0 (run El, short-dashed line), 2.5 
(run E2, long-dashed line) and 5 (run E3, dot-dashed line). The 
initial density profile (solid line) is plotted for reference. All pro- 
files have been normalised to po, the initial core density. Best fit 
values for po and r^ are listed in Table pi 



Run no. 


^impact 


rcoro (kpc) 


Pc/PO 


Initial 


- 


108 


1 


El 





125 


0.85 


E2 


2.5 


135 


0.65 


E3 


5 


157 


0.46 



Table 2. Parameters for the best fits to Equation hi for the density 
profiles plotted in Figure H. 



creasing the core density, with a trend towards a lower core 
density for the mergers with a non-zero impact parameter, 
with the core gas receiving additional rotational support. 
The enlargement of the core is a result of an increase in 
the entropy of the core gas during the merger, which can 
be seen in Figure M. Here, we plot the mean value of the 
entropy-related parameter 

Sr = U/P-^' (8) 

of the 250 particles initially closest to the the core of one of 
the clusters during two head-on mergers, one in which the 
clusters have equal masses (run El) and one in which the 
second cluster has one eighth of the mass of the first (run 
Ul). The cluster cores start to interact at i = 0, and the 
entropy is normalised so that it is equal to 1 at t = — 6Gyr. 
Initially the entropy decreases slowly as a result of radiative 
cooling but jumps sharply at t = due to the strong shock 
formed as the cores collide. The equal mass merger then 
experiences a second, larger, jump in entropy at i ~ 6Gyr 
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Figure 6. Evolution of the mean entropy of the core gas (defined 
in the text) during two head-on mergers, one between two equal 
mass clusters (run El, solid line) and the other between two clus- 
ters with masses in the ratio M1/M2 = 8 (run Ul, dashed line). 
The entropy has been normalised relative to the mean entropy at 
t = -6Gyr. 



Figure 7. Evolution of the boloinetric X-ray luminosity, defined 
by Equation M, during the merger of two equal mass systems. Four 
simulations are shown, with impact parameters 6 = (run El, 
solid line), 2.5 (run E2, dotted line), 5 (run E3, short-dashed line) 
and 10. (run E4-, long-dashed line). Times have been set so that 
t = corresponds to the peak of the boloinetric X-ray luminosity 
in the head-on merger. 



as the dark matter cores of the subclusters turn around and 
recoUapse (see PTC94 for a detailed discussion of this effect, 
which is absent in their gas-only simulations), while the un- 
equal mass merger experiences only one shock. The increase 
in core entropy at late times is due to the dissipation of 
oscillations as the core settles down after the collision. 



3.3 X-ray Luminosity, Temperature, and the 
Lx — 7x Relationship 

X-ray observatories allow precise measurement of the bolo- 
inetric X-ray luminosities (Lboi) and emission- weighted tem- 
peratures (Tew) of clusters of galaxies, which are observed 
to be correlated with approximately Lboi tx T^„ (Edge 
& Stewart 1991; David et al. 1993). This relationship is 
in conflict with theoretical models that assume clusters 
form through a self-similar gravitational collapse, which pre- 
dict M oc Tx^^ (Horner, Mushotzky & Scharf 1999) and 
Lhoi tx rx^(l + 2/)^/^ (Scharf & Mushotzky 1997) where Zf 
is the redshift at which the cluster forms, indicating that 
the ICM does not share the self-similarity seen in dark- 
matter only models. More recent work shows signs of a con- 
vergence between theory and observation, as correcting the 
observations for the effects of cooling flows flattens the re- 
lationship at high temperatures (Markevitch 1998; Allen & 
Fabian 1998), while more realistic numerical simulations in- 
cluding radiative cooling steepen the theoretical predictions 
(Muanwong et al. 2001). In principle, the dependence on Zf 
means that the observed Lboi — Tow relationship contains 
information about the evolution of clusters. However, while 



the observed relationship contains significant scatter (Allen 
& Fabian 1998), some of this is due to physical processes in 
the ICM such as radiative cooling (Pearce et al. 2000; Muan- 
wong et al. 2001) and possibly preheating (e.g. Lloyd-Davies, 
Ponman & Cannon 2000) which are not yet well understood. 
In addition merger shocks will also alter the observed lumi- 
nosity and temperature of clusters, and we examine these 
effects here. Throughout this section, we calculate the bolo- 
metric luminosity of the merging cluster as 



Lboi = / 



rriipi 



:MT„Z), 



(9) 



where the subscript i denotes the sum over all gas particles, 
which have temperatures Ti, densities pi and masses rrn. 
We assume a mean molecular mass /xmn = 10~^*g, and 
an emissivity A{Ti,Z) tabulated by Sutherland & Dopita 
(1993) with a metallicity Z — 0.3^©, where Zq is the solar 
value. The emission-weighted temperature is calculated as 



J cw — 



'^.■mipik{Ti, Z)Ti 



(10) 



^.miPiA(Ti,Z) • 

Figure W plots the evolution of the bolometric X-ray 
luminosity during simulations El — E4, which follow the 
merger of two equal-mass clusters at a range of impact pa- 
rameters. Times are scaled so that t = corresponds to the 
maximum brightening during the head-on merger. The peak 
luminosity occurs when the cores of the two clusters inter- 
act, and, in the case of head-on and nearly head-on impacts, 
a large increase in luminosity is apparent, with Lboi increas- 
ing to more than 5 times its precoUision value over a period 
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Figure 8. The evolution of the bolometric X-ray luminosity dur- 
ing two head-on mergers of equal-mass systems, carried out both 
with radiative cooling turned on (solid line) and with it turned 
off (dotted line). Times are scaled as in Figure M 



of roughly one sound crossing time. Prior to this time no 
increase in the total X-ray luminosity is apparent, despite 
there being clear evidence for a merger being in progress 
in temperature maps of the cluster (e.g. Figure Q, panel 3). 
After i = 0, the core of the cluster undergoes a period of 
expansion driven by the dark matter, leading to a decrease 
in the core density and a large (~ 10 x) decrease in X-ray 
luminosity, before experiencing a small secondary bounce 
(which occurs at f ~ 3 for the head-on merger, and slightly 
later for the simulations with a non-zero impact parameter) 
and settling into hydrostatic equilibrium. The final X-ray 
luminosity is lower than the sum of the emission from the 
two sub-clusters due to the increase in the entropy of the 
core gas during the merger. This conflicts with the results 
of Ricker (1998), who find an increase in the final luminosity 
after the merger, but this is likely to be an artifact of the lack 
of dark matter in those simulations, and results from Ricker 
& Sarazin (2001) appear to be close to those presented here. 
The total X-ray luminosity of the merger remnant is greatest 
in the head-on merger, with a trend towards a lower final lu- 
minosity with increasing impact parameter in keeping with 
the incr ease in core entropy with impact parameter seen in 
Section 



3.2 



The effect of including radiative cooling can be seen in 
Figure H, which plots the evolution of the bolometric X- 
ray luminosity in run El against an identical merger in 
which radiative cooling is not included. Radiative cooling 
increases the peak bolometric luminosity during the merger 
by roughly 20%, and leads to the gas in the core of the 
post-merger cluster having slightly lower entropy and conse- 
quently being more luminous than the gas in the simulation 
without radiative cooling. 



> 

CD 







Figure 9. Evolution of the X-ray emission- weighted temperature, 
defined by Equation hol during the merger of two equal mass 
systems. Four simulations are shown, with impact parameters b = 
(run El, soUd line), 2.5 (run E2, dotted fine), 5 (run E3, short- 
dashed line) and 10. (run E4, long-dashed line). Times are scaled 
as in Figure M. 

Figure plots the emission- weighted temperature of the 
cluster during the same set of simulations. Like the X-ray 
luminosity, the emission- weighted temperature increases sig- 
nificantly during the merger with the effect once again being 
particularly marked for mergers with a small impact param- 
eter, for which Tew increases by a factor of approximately 
3.5. Adiabatic cooling during the core bounce is clearly visi- 
ble, with the cluster then settling down at a higher temper- 
ature as would be expected given the higher virial mass of 
the remnant. 

The large changes in the X-ray luminosity and emission- 
weighted temperature of the cluster during the merger has 
implications for the Lboi — Tew and M — Tew relationships. 
Figure hd plots the evolution of the cluster in the Lboi — 
Tew plane during two equal mass mergers, one head-on and 
the other off-centre with 6 = 5. The dashed line in Figure uu 
represents a power-law fit to the Lboi — Tiw relationship in 
the simulations of Muanwong et al. (2001). The net move- 
ment on the Tew — iboi plane is not that great, given that the 
initial luminosity plotted in Figure flG is twice the value for 
an individual subcluster. The cluster becomes hotter due 
to the increased virial mass, but the injection of entropy 
into the core limits any increase in the luminosity of the 
system. However there are large movements when the sub- 
cluster cores merge, with the system initially following a 
track with L oc T'^ as the core gas is compressed adiabati- 
cally with the track steepening as the gas shocks at i = 0, 
and during this time the cluster appears to be much hot- 
ter and more luminous than its pre-merger state. The clus- 
ter subsequently returns to a constant entropy track as the 
core expands after the merger, with both the luminosity and 
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Figure 10. The evolution of the merging system in the Tew — 
Lboi plane during two equal-mass mergers, one head-on (run El, 
solid line) and the other with fe = 5 (run E3, dotted line). The 
dashed line represents a power-law fit to the simulations of Muan- 
wong et al. (2001) with Lbol <^ ^"^cw The high initial luminosity 
is a result of plotting the sum of the X-ray emission from the two 
subclusters. 



Figure 11. The evolution of the merging system in the Tow — 
Lbol plane during two mergers between clusters with masses in 
the ratio M1/M2 = 8 , one head-on (run Ul, solid line) and the 
other with 6 = 5 (run U3, dotted line). The dashed line represents 
a power-law fit to the simulations of Muanwong et al. (2001) with 
Lbol oc feTpw- 



temperature dropping below their pre-merger values before 
the core recoUapses. The implications for the M — Tow re- 
lation are also significant. The simulations of Muanwong et 
al. (2001) find that M oc {kT)^'"^ , and so a mass determina- 
tion for the cluster based upon this relationship would vary 
by a factor of ~ 9 depending on when the observation was 
taken. The strong variations in temperature and luminos- 
ity are relatively short lived, taking place on a timescale of 
around a sound crossing time, but for large clusters this can 
still represent lO' years or more. 

While equal-mass mergers will cause large fluctuations 
in the X-ray luminosity and the emission-weighted temper- 
ature of a cluster of galaxies, they are rare events. Unequal 
mass mergers will be much more common, and the increase 
in Lbol and Tow during these events may prove to be a sig- 
nificant source of scatter in the Lboi — Tow relationship. Fig- 
ure [nl plots the evolution of the cluster in the Tow — iboi 
plane during the merger of two systems with masses in the 
ratio 8:1. The maximum increase in luminosity is smaller 
than that seen in the equal-mass mergers, as would be ex- 
pected, but can still double during a head-on event. The 
lower post-merger luminosity indicates that there is still an 
injection of entropy into the core of the larger cluster, while 
the emission-weighted temperature of the cluster increases 
by as much as 2keV during the head-on merger. 



correspondingly long cooling times. However in ~70% — 90% 
of observed clusters the central gas has a radiative cooling 
time less than a Hubble time (Edge, Stewart & Fabian 1992; 
White, Jones & Forman 1997). This short cooling time leads 
to a slow inflow of gas to maintain pressure support known 
as cooling flow (Fabian 1994), in which as much as 1000 
M0yr"^ can cool out of the ICM (Allen et al. 1996). There 
is a significant anticorrelation between substructure in clus- 
ters and the presence of cooling flows (Buote & Tsai 1996) 
which are almost never associated with very irregular clus- 
ters (Edge et al. 1992), but the widespread nature of cooling 
flows suggests that they cannot be easily disrupted by minor 
mergers, which occur relatively frequently. 

The short cooling times and high X-ray luminosities in 
the cores of cooling flow clusters imply high gas densities. 
To assess the effects of this, we have carried out a second 
series of simulations of mergers of clusters containing cooling 
flows. These simulat ions are set up in a similar way to that 
described in Section 



2.1 



but with the core radius Vc reduced 
to 40kpc to give a higher core density. The outer density 
cutoff remains 1.6Mpc from the centre of the cluster, and 
so now R = 40rc. This density profile is no longer static 
as the cooling time is now comparable to the timescale of 
our simulations, and so we allow the clusters to evolve in 
isolation until an approximately steady state is reached (i.e. 
the density profile in the core has stopped evolving) before 
bringing them together prior to merging. The density profile 
of these 'cooling flow' clusters after this initial period of 
relaxation is shown in Figure 121 with the core density being 
with relatively low central densities (po = 10~"^cm~'^) and approximately 0.02cm~^ and the central cooling time being 



3.4 Cooling Flows 

The simulations in Section 



3.3| are concerned with clusters 
3 
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Figure 12. Radial density profile for our 'cooling flow' cluster 
(dashed line), with the initial density profile plotted for compar- 
ison (solid line). Both profiles have been normalised to po, the 
initial core density in the low-density cluster, and Tc, the core 
radius in the low-density cluster. 



around 2.5 x 10^ years. Prior to the merger the cooling flow 
is depositing mass at a steady rate of just over 8OM0yr~^. 
The evolution of the bolometric X-ray luminosity with 
time, shown in Figure |lj, is noticeably different from the low 
core density mergers examined in the previous section. The 
two X-ray brightness curves have been normalized to their 
luminosity at i = — 6Gyr for comparison, as the cooling-flow 
cluster otherwise has a much higher luminosity due to its 
higher core density. Unlike the low core density cluster, the 
cooling flow cluster brightens only briefly, although the peak 
luminosity is only slightly lower than in the low density sim- 
ulation. In addition, the second core bounce happens much 
more rapidly, and the cluster dims by a factor of more than 
30 between bounces. These effects are a result of the differ- 
ent density profiles of the two clusters. The X-ray luminosity 
of the two clusters is dominated by emission from the dens- 
est core gas, and will increase when that gas is compressed. 
The strongly-peaked density profile of the cooling-flow clus- 
ter will therefore increase its luminosity when the very centre 
of the two clusters interact, whereas the flat, constant den- 
sity core of the low-density cluster will experience a much 
more prolonged increase in luminosity. In addition, the more 
centrally concentrated mass distribution in the cooling flow 
cluster will relax more efficiently, leading to a shorter time 
between core bounces (the same effect can be seen in Pearce 
et al. 1993, who examine mergers between coUisionless sys- 
tems having Hubble density profiles with s=2, 3 and 4). 
The path of the cooling flow cluster in the Tow — iboi plane 
is similar to that seen in the previous section, although there 
is continuing evolution at late times as the dense core gas 
radiates away the energy gained in the merger. 



Figure 13. Evolution of the bolometric X-ray luminosity, de- 
fined by Equation bl during the head-on merger of two equal 
mass systems. The solid line represents the merger of two systems 
with core densities of 10~^cm~^ (run El"), while the dashed line 
represents the merger of two 'cooling flow' clusters with central 
densities of 10~^cm~'' (run CI). The luminosities have been nor- 
malized to their values at t = — 6Gyr for comparison, as the high 
density of the gas in the core of the cooling flow cluster mak es it 
much more luminous than the clusters examined in Section p.3[ 
Times are scaled as in Figure M. 



The unequal-mass merger in which the major clus- 
ter has a high density core show little evolution in the 
Tow — iboi plane. The head-on merger (run C^) shows a very 
brief increase in luminosity of roughly a factor of two and 
an increase in the emission-weighted temperature of around 
1.5keV, while an off-centre merger with & = 2.5 (run C5) 
shows only a small (^^ 10%) increase in luminosity and a 
barely-noticeable increase in the emission-weighted temper- 
ature (~ 0.2keV). The implication is that these unequal- 
mass mergers have only a small effect on the core gas of cool- 
ing flow clusters, especially when the merger is off-centre, 
and this can be clearly seen in Figure [l3, in which the evo- 
lution of the entropy of the gas in the core of the cooling 
flow cluster is plotted. In the head-on equal-mass merger 
the core gas receives a large increase in entropy, suggest- 
ing that the whole core has been disrupted, and the cooling 
time in the core jumps to around 10^° years. The actual 
mass deposition rate, measured from the number of gas par- 
ticles cooling to lO^K each timestep, drops to zero as the 
cores of the two clusters collide, increasing slowly after the 
merger, although even after lOGyr the mass deposition rate 
is still less than 2OM0yr~^. However, the apparent mass de- 
position rate derived from the luminosity and temperature 
within the cooling radius 
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Figure 14. The evolution of the mean entropy of the core gas 
(defined as in Figure H) during three mergers between 'cooling 
flow' clusters, including a head-on merger between equal-mass 
systems (run CI, solid line) and two mergers between two clusters 
with masses in the ratio M1/M2 = 8, firstly head-on (run C2, 
dotted line) and then with an impact parameter b = 2.5 (run C3, 
dashed line). Each run has been normalised relative to the mean 
entropy at t = — 6Gyr. 



(Fabian 1994) increases briefly during the merger and does 
not provide an accurate estimate of the mass cooling out 
of the ICM until the cluster has returned to hydrostatic 
equilibrium. 

The head-on unequal-mass merger also shows a jump 
in entropy, and although the cooling time of the core gas 
doubles to around 4 x 10^ years the energy injected into 
the core rapidly starts to be radiated away. The mass de- 
position rate again drops sharply during the merger, but 
recovers to around 55MQyr~^ within 4-5 Gyr. In contrast 
to the head-on mergers, the unequal-mass off-centre merger 
does not display any signs of the core gas being shocked, in- 
dicating that the bow-shock around the infalling subclump 
cannot have penetrated all the way into the core of the clus- 
ter. However the entropy of the core gas is decreasing more 
slowly after the merger, reflecting a slight increase in the 
coohng time to around 3 x 10^ years accompanied by a drop 
of around 30% in the mass deposition rate. The increase in 
the central coohng time is a result of ram pressure from the 
infalling subcluster displacing gas from the core of the clus- 
ter (Fabian & Daines 1991; Gomez et al. 2000). Ram pressure 
also accounts for the bulk of the disruption to the head-on, 
unequal-mass merger, as the increase in core entropy alone 
is insufficient to account for the doubling of the cooling time 
in the cluster core. 

The robustness of high-density cores to the effects of 
what should be a relatively common merger event is of in- 
terest to the survival of cooling ffows in clusters of galaxies. 
Our simulations suggest that a major merger between clus- 



ters of equal mass will completely disrupt any cooling ffow 
activity, but cooling flows, once established, may prove re- 
silient to mergers with lower-mass subclumps, especially if 
they are off-centre. While the cooling time in the low core 
density simulations remains greater than a Hubble time dur- 
ing almost all of the merger process, dipping below 10^" 
years only during the maximum compression of the cores, 
the cooling time in the high-core density clusters remains 
below a Hubble time throughout the process. In the case of 
the equal-mass merger the increase in core entropy and cool- 
ing time is large, and is enough to both disrupt the cooling 
flow and to prevent it from restarting soon after the merger, 
while the increase in the cooling times in the unequal-mass 
mergers are much less, and would indicate that the cooling 
flows would either experience minimal disruption or would 
rapidly be reestablished. 



3.5 Mixing 

Observations of the abundance of metals in the ICM sug- 
gest that negative radial metallicity gradients are common 
in clusters of galaxies (Irwin & Bregman 2001 and references 
therein). Significant gradients are most common in clusters 
containing cooling fiows, while clusters with no cooling fiow 
also tend to have a flat metallicity gradient, although this 
is not universal. The metallicity gradient is thought to be 
a result of either galactic winds or the stripping of metal- 
enriched gas from galaxies, and Irwin & Bregman (2001) 
suggest that the non-cooling-fiow clusters have experienced 
a merger that both disrupted the cooling flow and mixed 
the ICM, thus erasing the metallicity gradient. 

We examine the degree of mixing in the merger remnant 
in our simulations by smoothing the gas particle distribution 
on to a grid, and measuring a quantity A4 in each cell, where 



M 



1 



ni 



■ 712 



Wi +n2 



(12) 



and ni and n2 are the number of particles in the cell which 
were in cluster 1 and 2 respectively at the start of the simu- 
lation. The cell size is set so that each cell will, on average, 
contain at least ten particles, with the cells in the core con- 
taining more than this number, although this is not always 
true in the low-density outer regions of the cluster which are 
represented by relatively few particles. In a well-mixed cell 
ni and n2 will be equal and so A^ ^ 1, while in a cell in 
which no mixing has occurred A4 ~* 0. 

The degree of mixing 5Gyr after a head-on, equal-mass 
merger is shown in Figure [13. Very little mixing of the ICM 
has taken place, with the mixing in the centre of the cluster 
only appearing at around t = 4Gyr, apparently driven by 
the infall of material back along the merger axis rather than 
by the merger itself. Much more mixing is visible in Fig- 
ure Il6l which shows the state of an off-centre merger with 
& = 5. This is also plotted at t = 5Gyr, although in this case 
the cluster is still settling back to a steady state, as can be 
seen in the distortion of the X-ray isophotes. The angular 
momentum in the system means that the core of the merger 
remnant is well mixed, although there is still little mixing 
of the ICM at r > 6. 

Finally, Figure hj shows the degree of mixing during 
an off-centre merger in which the main cluster has a mass 
eight times that of the infalling subcluster. In this case, low- 
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Figure 15. Mixing 5Gyr after a head-on merger between two 
equal mass clusters (run El ) . Grayscales show the degree of mix- 
ing (defined in the text), contours represent logarithmic only X- 
ray isophotes and arrows represent the velocity field. 



Figure 17. Mixing 5Gyr after a off-centre merger with impact 
parameter b = 5 between two clusters with masses in the ration 
Mi/M2_= 8 (run US). Grayscales, contours and arrows are as in 
Figure ha. 




Figure 16. Mixing 5Gyr after a off-centre merger with impact 
parameter fe = 5 between two equal mass clusters (run E3). 
Grayscales, contours and arrows are as in Figure llSl 



density gas is rapidly stripped from the subcluster and mix- 
ing occurs throughout the outer layers of the main cluster. 
However, the subcluster core is less efficiently stripped, and 
can be seen spiraling all the way into the core of the main 
cluster. There is a clear difference between the final states 
of the equal mass and unequal mass mergers. Very little 
mixing occured in the outer layers of the cluster during the 
equal mass mergers, although the off-centre mergers mixed 
the cores efficiently, while the unequal mass merger caused 
very little mixing in the core but lead to significant mixing 
in the outer regions of the cluster. 

While none of these mergers mixed the ICM globally, 
our simulations suggest that mergers of the magnitude nec- 
essary to disrupt a cooling flow can efficiently mix the cores 
of clusters on the scales accessible to X-ray observations (in 
general < 50% of the virial radius) if the merger is off-centre. 
In contrast, the unequal mass merger was unable to mix the 
core gas efficiently, but may cause locally increased abun- 
dances such as that seen by Arnaud et al. (1994) in the 
Perseus cluster and are consistent with the presence of a 
cooling fiow, which will not necessarily be disrupted during 
a minor merger. 



3.6 Numerical Resolution 

Quantifying the effects of limited resolution is always im- 
portant in numerical simulations. To assess the impact that 
resolution has on the results presented here, we have car- 
ried out two further simulations of head-on mergers between 
equal mass systems in which the total number of particles 
have been reduced by a factor of 4 (run Rl ) and 16 (run 
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Run no. 



JVtot mdm(MQ) mgas(MQ) 



El 131072 1.25 X 10^° 1.13 x 10^ 

Rl 32768 5 x 10^° 4.5 x 10^ 

R2 8192 2 X 10^1 1.8 x 10^" 

Table 3. Particle numbers and masses for the simulations dis- 



cussed in Section 3.6 
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Figure 18. Evolution of the merging system on the Tew — 
Lbol plane during during three head-on mergers between equal 
mass systems. The solid line represents a simulation in which each 
cluster has 32768 particles each of gas and dark matter (run El), 
the dotted line a simulation with 8192 particles of each type (run 
Rl ) and the dashed line a simulation with 2048 particles of each 
type (run R2). The total mass of the systems is kept constant by 
increasing the particle masses in the lower resolution runs. 



R2), giving a total of 32768 and 8192 particles respectively, 
with the cluster properties kept constant by increasing the 
mass of the remaining particles. Details of the simulation 
parameters are shown in Table pi 

The evolution of the clusters in the Tew — iboi plane dur- 
ing these simulations is shown in Figure hq. The convergence 
between runs is good, with the evolution of the system being 
similar in all three simulations. The main effect of lowering 
the resolution is a systematic decrease in the luminosity of 
the merger remnant. In our low-resolution runs less energy 
is deposited in the halo gas than in the full runs and the 
core gas is shocked to a greater extent, leading to a broader, 
higher-entropy core in the final cluster, a result also seen in 
the resolution testing carried out in PTC94. The differences 
in the emission-weighted temperature are small. 

If the particle number is reduced further then the agree- 
ment between simulations worsens, and we would agree with 
the conclusions of PTC94 that 1000 gas particles per clus- 
ter represent a minimum for modelling the bulk properties 



of a merging system. However, more than 10^ gas particles 
per cluster are required to produce well-resolved tempera- 
ture and X-ray surface brightness maps such as those shown 
in Figures H and H, in agreement with Steinmetz & Miiller 
(1993). 

We conclude that the effects of resolution are of limited 
significance in our simulations. The evolution of the cluster 
in the Tew — Lhoi plane has converged to a consistent result, 
differing only in detail as the resolution is increased. Equally, 
we have sufficient resolution to follow the propagation of 
shocks in the ICM acceptably, although here our resolution 
is closer to the minimum required for the problem. 



4 DISCUSSION 

It is clear from the results presented here that mergers can 
alter the global properties of clusters of galaxies. During a 
merger, hydrodynamic shocks dissipate much of the kinetic 
energy of the impact into the ICM, leading to a departure 
from hydrostatic equilibrium that can last for several Gyr. 
Mergers lead to many potentially observable effects, includ- 
ing strong temperature gradients across shocks, distortion 
in the X-ray isophotes, multiple peaks in the X-ray surface 
brightness and high- velocity bulk fiows in the ICM. Despite 
this, identifying merging clusters is not necessarily easy, as 
both projection effects and the limits of X-ray observations 
can serve to mask many of the signs of an ongoing merger. 
While X-ray substructure is clearly visible both before and 
after a merger when viewed perpendicular to the merger 
axis (as can be seen in Figures 3 and 0), the X-ray isophotes 
appear spherically symmetric throughout the merger when 
viewed along the merger axis. Equally, the shocks so appar- 
ent in temperature maps when viewed perpendicular to the 
merger are not obvious when viewed along the merger axis, 
with the cluster just appearing to have a steep temperature 
gradient as the cores of the two subclusters merge. While 
this departure from isothermality is itself a good indicator 
of a merger, it requires that the temperature profile is both 
sufficiently well resolved for it to be apparent and is not 
unduly affected by the effects of the merger, as the X-ray 
deprojection method (e.g. Arnaud 1988) assumes that the 
cluster is spherically-symmetric and in hydrostatic equilib- 
rium. Temperature maps, derived directly from the X-ray 
hardness ratio, will therefore be a better indicator of an on- 
going merger than deprojected temperature profiles. If the 
cluster is poorly resolved, it will simply appear to be well 
relaxed, hot and luminous, although evidence for a merger 
will still be seen in the distribution of the peculiar velocities 
of the galaxies in the two clusters, which can be expected to 
be significantly bimodal during a merger. At intermediate 
angles, projection effects will serve to smooth distortions in 
the X-ray isophotes and reduce the visibility of the merger in 
the X-ray waveband, although the effects are still apparent. 
The merging process has major implications for statis- 
tical studies of clusters. As we have shown in Section |3.3| , 
even relatively minor (M1/M2 — 8) mergers can introduce 
a significant scatter in the observed Lboi — Tow relationship, 
while estimates of the mass of clusters based on the M — Tow 
relationship (Horner et al. 1999) will suffer large uncertain- 
ties during and immediately after a merger. Mass estimates 
based on assumptions of hydrostatic equilibrium (e.g. Fabri- 
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cant, Lecar & Gorenstein 1980) will also be strongly affected 
as the departures from equilibrium can last for several Gyr 
after the merger, during which time many of the readily- 
observed signs of the merger have faded, and Roettiger et 
al. (1996) estimate that errors in the mass estimate may 
reach 50% in the 2Gyr following a merger. Edge et al. (1992) 
suggest that clusters will typically undergo a merger every 
2-4Gyr, similar to the timescale for our clusters to return 
to hydrostatic equilibrium, implying that truly relaxed clus- 
ters may be uncommon. The large increase in X-ray lumi- 
nosity during a merger will also introduce strong selection 
effects into cluster surveys, and the clusters found at high 
redshift may not be a representative sample. The increases 
in X-ray luminosity and temperature during a major merger 
will strongly influence the statistics of the hottest and most 
luminous clusters of galaxies. This has been conflrmed by 
Chandra temperature maps of two of the hottest clusters 
known, IE 0657-56 (Tucker et al. 1998) and A2163 (David 
et al. 1993; Markevitch & Vikhhnin 2001), both of which 
show strong merger shocks. In addition, unusual clusters 
like A851, which is cooler and less luminous than it's rich- 
ness would suggest, may be undergoing the post- merger core 
expansion (Schindler & Wambsganss 1996). 

The situation is more complex when cooling flow clus- 
ters are considered. We find that major mergers will almost 
certainly disrupt cooling flows to the extent that they can- 
not restart within a Hubble time, consistent with the lack 
of cooling flows in irregular clusters (Edge et al. 1992), but 
our unequal mass mergers have less of an effect on the cool- 
ing flow. A head-on merger still disturbed the cooling flow, 
although it restarted within a relatively short period, while 
an off-centre merger had little impact on the flow. The sur- 
vival of the cooling flows during minor mergers is supported 
by a growing body of observations of cooling flow clusters 
that also show evidence for a late-stage merger (e.g. A2142; 
Markevitch et al. 2000) , as well as merging clusters with 
small cooling flow cores (e.g. A2065; Markevitch et al. 1999) 
which suggest that the cooling flow has been disturbed but 
has survived, at least until the present day. 

Finally, we examined the degree to which mergers can 
mix the ICM. While the ICM is never well-mixed globally, 
our simulations suggest that major off-centre mergers can 
efficiently mix the core gas, although there was very little 
mixing during a head-on merger. Equal mass mergers will 
also disrupt any cooling flow in the cluster, and may explain 
why metallicity gradients are rare in clusters without a cool- 
ing flow (Irwin & Bregman 2001). The unequal mass merger 
introduced a high degree of mixing in the outer layers of 
the main cluster as the low-density gas in the subclump is 
rapidly stripped away, but had little effect on the core gas. 

A feature of unequal-mass mergers that we have not 
touched on in this paper is the survival of the cores of sub- 
clusters during mergers. Sharp surface brightness disconti- 
nuities have been observed by Chandra in merging clusters 
of galaxies such as Abell 2142 (Markevitch et al. 2000) and 
Abell 3667 (Vikhlinin et al. 2001), and it has been suggested 
that these features, termed "cold fronts", are low-entropy 
gas from the core of a subcluster that had so far survived 
the merging process (Markevitch et al. 2000). We find that 
subcluster cores can survive off-centre mergers, and result 
in features very similar to those observed by Chandra ; this 
will be examined further in a subsequent paper. 



5 CONCLUSIONS 

We have presented results from high- resolution AP^M+SPH 
simulations of merging clusters of galaxies. Our principle 
conclusions are; 

(i) Mergers lead to many potentially observable effects, 
including strong temperature gradients across shocks, dis- 
tortion in the X-ray isophotes, multiple peaks in the X-ray 
surface brightness and high- velocity bulk flows in the ICM. 

(ii) An increase in the entropy of the core gas during the 
merger leads to a broader, lower-density core in the post- 
merger cluster. In addition, off-centre mergers can give the 
core additional rotational support. 

(iii) The compression and shocking of the core gas during 
a merger can lead to large increases in the bolometric X- 
ray luminosity and the emission-weighted temperature of 
the cluster. This will have a significant impact on statistical 
studies of galaxy clusters. Immediately after the cores of 
the subclusters interact, the core of the cluster undergoes a 
period of expansion driven by the coUisionless dark matter, 
cooling adiabatically and decreasing in luminosity by more 
than an order of magnitude, before the cores turn around 
and recoUapse. 

(iv) Radiative cooling leads to an increase of ~ 20% in the 
X-ray luminosity during and after a merger when compared 
to simulations without cooling, even if central densities are 
relatively low (po = lO^^cm"^) and the cooling times in 
the pre- and post-merger objects are greater than a Hubble 
time. 

(v) In a major merger the cluster can be out of hydro- 
static equilibrium for several Gyr, during which time the 
clear observational signs of a merger will become less appar- 
ent. Brief departures from hydrostatic equilibrium were also 
seen in simulations of unequal mass mergers. 

(vi) Cooling flows will be completely disrupted during 
major mergers, and are unlikely to be able to restart within 
a Hubble time. However, minor mergers had less of an effect 
on a cooling flow, which was either unaffected by the merger 
or was able to restart rapidly. Ram pressure is more effec- 
tive at disturbing the cooling flow in minor mergers than 
the merger shock. During a merger the mass deposition rate 
inferred from observations is a poor estimate of the actual 
mass cooling out of the ICM. 

(vii) Major off-centre mergers effectively mixed the gas 
in clusters within a radius roughly equal to b, the impact 
parameter. Head on mergers and unequal-mass mergers had 
little effect. None of the mergers examined here effectively 
mixing the ICM globally. 
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